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Several large-scale gravitational wave (GW) interferometers have achieved long term 
operation at design sensitivity. Questions arise on how to best combine all available 
data from detectors of different sensitivities for detection, consistency check or veto, 
localization and waveform extraction. We show that these problems can be formulated 
using the singular value decomposition (SVolHmethod. We present techniques based on 
the SVD method for (1) detection statistic, (2) stable solutions to waveforms, (3) null- 
stream construction for an arbitrary number of detectors, and (4) source localization for 
GWs of unknown waveforms. 

Keywords: Gravitational wave, data analysis, singular value decomposition 

1. Introduction 

Several large-scale gravitational wave (GW) interferometers have achieved long term 
operation at design sensitivity. Different detectors can have very different noise levels 
and different frequency bandwidth. The directional sensitivity of different detectors 
can also be very different. For instance, the most sensitive GW detectors in the 
US, namely the LIGO detector at Livingston, Louisiana, and the two co-located 
detectors at Hanford, Washington (abbreviated as LI and III/II2), are designed to 
be nearly aligned. Detectors in Europe (GEO600 in Germany, and VIRGO in Italy) 
and in Asia (TAMA in Japan) are incidentally nearly orthogonal in directional 
sensitivities to the LIGO detectors. The questions arise on how to best combine 
data from these detectors in GW data analysis. 

We present in this paper an application of the singular value decomposition 
(SVD) method^ to data analysis for a network of GW detectors. We show that the 
SVD method provides simple solutions to detection, waveform extraction, source 
localization, and signal-based vetoing. By means of SVD, the response matrix of 
the detector network can be decomposed into a product of two unitary matrices 
and a pseudo-diagonal matrix containing singular values. The unitary matrices can 
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be used to form linear combinations of data from all detectors that have one to one 
correspondence to linear combinations of the gravitational wave signal polarization 
components. Each newly formed data stream has a corresponding singular value 
representing the network's response to the new signal polarizations. Data streams 
with non-zero singular values represent the signal components while data with zero 
singular values (or zero multiplication factors) represent the null streams. The null 
streams have null response to GWs and can be used for localization of GW sourceJ^ 
and for identifying detector glitches from that of real gravitational wave events as 
proposed by Wen and Schutz (2005)'^ for ground-based GW detectors. The statis- 
tical uncertainty in estimating the GW waveforms from the data can be shown to 
be related to the inverse of the singular values. This can be used to reduce errors 
in signal extraction by enabling "bad" data with relatively small singular values to 
be discarded. 



2. Mathematical Framework 



The observed strain of an impinging GW by an interferometric detector is a linear 
combination of two wave polarizations 



hi{t) = fth+{t) + !l^hAt), 



(1) 



where t is time, /+ and are the detector's response (antenna beam pattern func- 
tions) to the plus and cross polarizations (/i+, ^x) of a GW wave. These antenna 
beam patterns depend on source sky directions, wave polarization angle, and de- 
tector orientation. The / subscript is a label for the /th detector, indicating the 
dependence of the observed quantities on detectors. 

Data from a GW detector can be written as the sum of the detector's response 
and noise, di{t + tu) = hi{t) + ni{t + tu), where, t = [0,T], T is the duration of 
the data used, tu is the wavefront arrival time at the /-th detector relative to the 
reference detector (labeled as detector 1). The arrival time delay can be calculated 
as Tu = h ■ rii/c, where ru is location vector of the /-th detector relative to 
the reference detector, c is the speed of light, fi is the unit vector along the wave 
propagation direction. 

For a given source direction, the time-delay corrected data from a network of 
Nd GW detectors can be written in the frequency domain as. 



dk = Akh 



feiik 



{k = l,...,Nk) 



(2) 



where, k is the label for frequency bins, Ak is the response matrix of the detector 
network at each frequency. 



Ik = 



/ dlk/cTlk \ 
d2k/(^2k 

KdNak/cTNak/ 



l-xk 



,nk 



/ nik/aik \ 
n2k/cr2k 

\nNik/<7NakJ 



(3) 
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and 



Ak 



(4) 



where cr?j. is the noise variance of the i-th detector at the fc-th frequency bin. The 
response of the detector network to GWs in all frequencies can be further written 
in terms of vectors and matrices as, 



where, 



d2 



h2 



VhNj 



d = Ah + n, 



n2 

,n 



(5) 



VnNk/ 



v 



(6) 



A of dimension N^Nk x 2Nk is therefore the response matrix of the detector network 
for all frequencies. 



3. Principle of Data Analysis Based on SVD 

The singular value decomposition of the network response matrix A can be written 



as. 



usV 



(si 

S2 

... 

... 



\0 



\ 


S2Nk 









(7) 



where singular values Si > S2 > ••• > S2Nk > 0. w and v are unitary matrices of 
dimensions N^Nk x NdNk and 2Nk x 2Nk respectively with uhi = I, v''v = I. Note 
that sf arc eigenvalues of A^A. Note also that the singular values of A are only a 
rearrangement of singular values of A^. 

We propose to construct new data streams by linearly recombining data from 
different detectors and construct new signal polarizations using the unitary matrices 
resulting from the SVD of the network response matrix A, 



d' = uM, h' = v'^h, 



n' = u^n. 



(8) 
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These newly constructed data streams and new signal polarizations are one-to-one 
related to each other via the pseudo-diagonal matrix containing singular values, 



/ 4 \ 
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(9) 



At the presence of a GW signal, there are at most two data streams at each fre- 
quency corresponding to non-zero singular values. If the data contain pure station- 
ary Gaussian noise only, the new data stream also follows Gaussian distribution 
and are statistically uncorrelated with each other. These data streams can therefore 
be treated effectively as one-detector data streams containing signals. All detection 
algorithms for a single detector can be applied to these "signal" data streams. 

The singular values of the response matrix A of the detector network are directly 
related to the statistical uncertainties in estimating waveforms from the data. The 
Fisher information matrix for estimating wave polarizations h is, 



/(h) 



Vo 







(10) 



where we define the variance of a complex number to be the sum of the variance 
of the real and complex parts. The lower-bound^ for the variance of the estimated 
GW signals due to statistical errors are therefore a linear combination of s~^. 
Singular values are therefore indicators for the accuracy of waveform estimation and 
can be used for regularization of the solutions to the waveforms. 

At the end, data streams corresponding to zero singular values or zero multipli- 
cation factors have null response to signals (Eq. [9]). There are at least {Nd — 2)Nk 
"null streams" due to the fact that there are only two polarizations for a plane GW 
based on Einstein's theory of General Relativity. Null streams can therefore be used 
to test^the consistency of the detected GW events and to veto against signal- like 
noise 



El 



3.1. Detection Statistic 

In this section, we discuss the construction of detection statistic for GWs of un- 
known waveforms using data from a network of detectors. We discuss two possible 
approaches, (1) a direct application of maximum likelihood raticP^ (MLR) principle 
by constructing detection statistic from adding powers of all data with non-zero 
singular values, (2) optimized detection statistic based on MLR together with as- 
sumptions on wave spectrum and principle of maximum signal-to- noise ratio (SNR). 



7, 2008 4:25 WSPC/INSTRUCTION FILE msastrodv2 



Detecting Gravitational Waves Using Detector Arrays 5 

SNR is a good indicator for detection efficiency in our case as, wfien adding sufficient 
large number of noise powers together, central limit theorem states that these SNRs, 
regardless of underlying noise distributions, follow roughly Gaussian distribution. 

The "standard" solution is based on a direct application of the MLR principle 
for unknown waveforms. When the signal polarizations /i+fe and hxk are assumed 
to be independent variables for all frequencies, the components of the new signal 
h' are also independent from one another. The standard MLR detection statistic 
can then be constructed by adding equally all powers of "signal" data-streams of 
non-zero singular values, 

Pf=EK:l'' (11) 

1=1 

where Np is the number of non-zero singular values. Note that singular values have 
already been ranked (Eq. [7]). In the absence of noise, jd^p = s^\h'\f. In the pres- 
ence of pure stationary Gaussian noise, 2|(i^p follows distribution with variance 
of 4. By discarding data of null response to GW signals, this detection statistic is 
optimal in SNRs than simply adding powers of all (noise- weighted) data together. 
However, in the presence of extremely small singular values (e.g, when two detec- 
tors are nearly aligned) or extremely weak signals, this detection statistic does not 
necessarily maximize the SNR. 

Singular values can be used to further optimize the detection statistic by making 
reasonable assumptions about the signal spectrum. For instance, we can omit signal 
powers of small singular values with hopes that the chance is very small for the new 
signal component h' to be very strong at frequencies and polarizations where the 
network is not sensitive. Under this assumption, the detection statistic can be, 

K 

pW^^ld'.p. (12) 

t=i 

That is, only < Np number of powers of reasonable large singular values are 
included. 

If we assume that the detection statistic is a linear combination of individual 
powers and that the signal power is the same in both frequency and new 

polarizations, maximization of the SNR leads to, 

K 

Pf =5]a.|(d').P, (13) 

i=l 

where, 

a, = ■ (14) 

If we assume flat power in signal polarizations only, summation and normalization 
of sf terms in Eq. [14] apply only to the two polarizations within each frequency. 
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We have also investigated the approach of adding powers within a window of 
the ranked singular values as a general approach to unknown waveforms, 

il+N' 

pf^ = E Ki'- (15) 

i—ii 

where the bandwidth of A^' is determined by maximizing the SNR, 

SNRj^, = max^, (P^v' /^)- (16) 

This is an ad hoc approach that works better for signal powers that are narrow- 
banded in frequencies. While this approach can be time-consuming, it performs 
pretty well in our test with the BH-BH merger signals. 

Early work on optimization of the detection statistic based on directional sensi- 
tivity of GW detectors can be found in Rajesh Nayak et al. (2003) for future spa. ce 
detector LISA^Iand in Klimenko et al. (2005, 2006) for ground-based detectors.^^^It 
has been demonstrated in these papers that further optimization from the standard 
MLR detection statistic can be obtained by constraining or eliminating contribu- 
tions of data corresponding to weaker network sensitivity. Other work can be found 
in Mohanty et al. (2006)P which is based on the SNR variability regulator and in 
Summerscales et al. (2006)P^ which is based on the maximum entropy principle. 

We show that the SVD method provides a simple mean to the construction of 
the detection statistic and its further optimization. The standard detection statistic 
based on the maximum likelihood ratio principle can be easily reproduced from 
the newly constructed data. Various optimization strategies can then be obtained 
based on general assumptions of the signals together with the singular values which 
represent the sensitivity of the detector network. We have also proposed a strategy 
by maximization of the SNRs from the summation of powers corresponding to a 
window of ranked singular values. A comparison o f th e performance of different 
detection statistic can be found in a follow-up paper.El 

3.2. Stable Solutions to Waveform Extraction 

The "standard" estimator for h' = [h[, according to the maximum likeli- 

hood ratio principle of "burst" GWs (defined as any GWs of unknown waveforms) 
is 

h\^s-^d\, i^l,...,Np, (17) 

where Np is the number of non-zero singular values. Note the components of h' 
include two polarization components for each frequencies. The "standard" solutions 
to wave polarizations can be extracted using the SVD components of the network 
response (Eq. [71 Eq. ^ 

h(") = vh'. (18) 

Note that the solutions of two signal polarization components at each frequency 
depends on quantities within that frequency only. Faster calculations can therefore 
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be carried out independently at each frequencies. In the hmit when Np = 2Nk, 
the solution given in Eq. [18] is equivalent to that written with the Moore-Penrose 
inverse, h(o) = (Al'A)-Mtd. Note that this is also the same as the effective one- 
detector strain for data from a network of GW detectors introduced by Flanagan 
and Hughes (1998).i2 

The "standard" solution h'^°) from Eq . [T8l however can possibly lead to unstable 
solutions in the sense that a small error in the data can lead to large amplified 
error in the solution. Fisher information matrix (Eq. I10|) indicates that the best 
possible statistical variance for the estimated is a linear combination of 1/sf. 
The extracted wave polarizations /iJJ can contain large errors if we include data 
corresponding to very small singular values. This situation occurs when the response 
matrix A is "ill-conditioned" with Si/si <C 1. The small singular values can result 
from machine truncation errors instead of zero values or from nearly degenerated 
solutions to the equations, e.g., in our case, when antenna beam patterns of two 
detectors are nearly aligned. 

Regularization is needed in order to have stable solutions to h^. One simple 
solution is to apply the "truncated singular value decomposition" (TSVD) method 
by omitting data with very small singular values in Eq. 1181 

Nt 

= Y,u,,h'^, (19) 

where iV^ < Np is the number of data included. The main problem is then the 
decision on where to start to truncate the data streams with small singular values. 
It depends on the accuracy requirement in waveform extraction, type of waveforms 
and type of constraints that can be put on the solutions. The fractional error due to 
truncation, defined as the ratio of the sum of error-square and the sum of detector- 
response square from all frequencies, is, 

Phil E,t^^|l(-th),p ^ ^ 

Truncation of terms with very small singular values can therefore retain the least 
square fit of the detector response to the data while greatly reduce the statistical 
errors when extracting individual signal polarizations. A recent work on introducing 
the Tikhonov regularization to waveform extraction of GW signals can be found 
in Rakhmanov (2006). Note that the Tikhonov regularization is equivalent to 
introducing a new parameter to filtering out data associated with small singular 
values. 

3.3. Null Stream Construction 

There are at least [Nd — 2)Nk data streams with zero singular values or zero mul- 
tiplication factors. These are data streams that have null response to signals. The 
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"standard" null streams can be written in terms of the SVD of A (Eq. [7]) as, 

^{u^d)^+N^, i = 1,2,.. NdNk ~ Np. (21) 

For stationary Gaussian noise, follows Gaussian distribution of zero mean and 
unity variance and are statistically independent with each other. Consistency check 
for a detected GW event therefore requires that all the null-streams are consistent 
with 'noise-only" at the source direction. Using null-streams as a tool for consis- 
tency check of GW events against signal-like glitches for ground-based GW de- 

o 

tectors was first proposed by Wen and Schutz (2005). Further investigations on 
consistency check using the null-streams and its relation to the null space of the 
network-response matrix A can be found in Chatterji et al. (2006). 121 

There are also "semi-" null streams corresponding to data streams with very 
small but non-zero singular values, 

Nw^ = iu^d),+N^^l, z = l,2,...,7Vp-iV;, (22) 

where is the starting indexes for 1 > 0. These semi-null streams exist when 
the equations are nearly degenerated. This happens, for example, for response of 
the two LIGO detectors of HI (H2) and LI that are designed to be nearly aligned. 
The semi-null stream can be also simply caused by combination of weak directional 
sensitivity and/or high noise level instrument of all detectors in the network. Con- 
sistency check of GW events, veto against noise, and localization of the source can 
be further improved by including both the "standard" null-streams of analytically 
zero singular values and also these semi-null streams. 

4=1 1=1 

The usage of semi-null streams depends on waveforms and therefore efficiency stud- 
ies should be carried out before-hand. A proposal of using the semi-null stream for 
veto against detector glitches and for source localization can be found in Wen and 
Schutz (2005)Pfor the two-detector network of Hl-Ll. 

4. Conclusion 

A new data analysis approach based on the singular value decomposition method 
has been proposed for the data analysis of GWs using a network of detectors. We 
show that singular values of the response matrix of the GW network directly encode 
the sensitivity of the GW network to signals and the uncertainties in waveform es- 
timation. The SVD method is particularly useful for constructing null data streams 
that have no or very little response to GW signals. We argue that the SVD method 
provide a simple recipe for data analysis of GWs for a network of detectors. Note 
that the SVD method is widely used in engineering signal processing for image 
processing, noise reduction and geophysical inversion problems. An application of 
the SVD method to detecting GWs from the extreme-mass-ratio-inspiral sources 
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using the space GW detector LISA can be found in Wen et al. (2006)!IS1 The SVD 
software package is also available in MATLAB. 

Strategies on construction of the detection statistic, stable solutions to wave- 
forms, and null streams based on the SVD method are discussed. We show that 
the detection statistic should be constructed from data streams of non-zero singular 
values. We discuss how detection efficiency can be improved from a direct applica- 
tion of MLR by incorporating our knowledge of the network's directional response 
to GWs and our assumptions on distribution of the signal power. 

We also give expressions of null-streams for arbitrary number of GW detectors 
using components from the SVD of the network-response matrix. The concept of the 
semi-null streams that are characterized by small singular values is also introduced. 
We argue that the exploration of the null-space by including semi-null streams can 
help improving the source localization and consistency check. Analytical study for 
the angular resolution of a network of GW detectors and their relations with the 
null space will be found in WeiPl (2007). We also show how a stable solution to 
the waveform can be constructed based on information from the singular values. 

We conclude that a GW event should be identified only when both of the fol- 
lowing two conditions are satisfied at the same sky direction, (1) high probability of 
detection for the optimal statistic constructed from "signal" streams with non-zero 
singular values (Eq. [51 sec. 13. ip and (2) high probability that the null streams cor- 
responding to zero singular values are consistent with noise (Eq.[2 Eq. [211 sec. l3.3p . 
Results of an extended investigation will be published elsewhere. 
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